MADPH-10-1566 



Lepton Number Violation from Colored States at the LHC 

Pavel Fileviez Perez, Tao Han, Sogee Spinner, and Maike K. Trenkel* 
University of Wisconsin-Madison, Department of Physics 
1150 University Avenue, Madison, WI 53706, USA 
(Dated: December 10, 2010) 

The possibility to search for lepton number violating signals at the Large Hadron Collider (LHC) 
in the colored seesaw scenario is investigated. In this context the fields that generate neutrino masses 
at the one-loop level are scalar and Majorana fermionic color-octets of SU(3)c- Due to the QCD 
strong interaction these states may be produced at the LHC with a favorable rate. We study the 
production mechanisms and decays relevant to search for lepton number violation signals in the 
channels with same-sign dileptons. In the simplest case when the two fermionic color-octets are 
degenerate in mass, one could use their decays to distinguish between the neutrino spectra. We find 
that for fermionic octets with mass up to about 1 TeV the number of same-sign dilepton events is 
larger than the standard model background indicating a promising signal for new physics. 



I. INTRODUCTION 



The LHC Era brings with it the hope of the discovery of the New Standard Model of particle physics. 
Searching for lepton and/or baryon number violation at the LHC might be an integral part of establishing 
this new theory. As it is well known, in the Standard Model (SM) the lepton and baryon numbers are 
accidental global symmetries at the classical level. Typically, beyond the SM physics scenarios introduce 
new interactions which can violate these symmetries. A concrete example are Majorana neutrinos, which 
require extra states and interactions to generate neutrino masses. These can give rise to new and interesting 
phenomenology relevant for the LHC. 

Recently, a simple mechanism for the generation of Majorana neutrino masses at one-loop level was 
proposed where the new fields, inside the loop, live in the adjoint representation of SU(3)c CD- In this 
context two types of new fields are required: a fermionic octet, p ~ (8, 1,0), and a scalar octet S ~ 
(8, 2, 1/2). It is important to stress that this is the simplest mechanism at one loop allowed by cosmology 
which does not require an extra symmetry. We refer to this model as the "colored seesaw mechanism". This 
is a novel and interesting possibility which one can use to test the origin of neutrino masses and look for 
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lepton number violating signals at the LHC since the seesaw fields can be easily produced via the strong 
interaction. 

In this article, we investigate these lepton number violating signals at the LHC in the context of the 
colored seesaw mechanism. We study the production mechanisms and decays of the colored states, focusing 
on same-sign dileptons in the final state as an indicator of lepton number violation. In the simplest case 
when the two fermionic color-octets are degenerate in mass, their decays allow us to distinguish between 
the neutrino spectra. We find that for fermionic octets with mass up to a few TeV the number of same-sign 
dilepton events is larger than the SM background. While our study is only on the level of cross section times 
branching ratio, it is still a promising result for the testability of this mechanism at the LHC. 

The paper is organized as follows: In Section II we set the stage for this paper by summarizing the 
potential for observing lepton number violation in the simplest extensions of the SM with Majorana neutrino 
masses. In Section III we describe in detail the colored seesaw mechanism and the constraints coming from 
neutrino physics and lepton flavor violation. Section IV is devoted to the discussion of the testability of the 
model at the LHC: we discuss the possible two- and three-body decays, and the production of the fermionic 
octets. We finish by summarizing our findings in Section V. 

II. LEPTON NUMBER VIOLATION AT THE LHC 

It is well known that neutrinos are either Dirac or Majorana fermions. Majorana neutrinos allow for a 
great variety of scenarios in which neutrino masses originate from integrating out heavy fields, the generic 
seesaw mechanism. A tantalizing possibility is that the new heavy fields could be produced at the LHC and, 
due to their Majorana nature, give rise to lepton number violating signals. 

At tree level, neutrino masses can be generated through the well-known Type I, Type II or Type III 
seesaw mechanisms: 

• Type 1 seesaw mechanism 10: The SM is extended by at least two SM singlets, v ~ (1, 1, 0). Once 
those singlets are integrated out the neutrino mass matrix reads as Ail = Y v M^ 1 Y^v 2 , where Y v 
is the Yukawa coupling between the SM leptonic doublet and the right-handed neutrinos, v = 246 
GeV is the vacuum expectation value of the SM Higgs boson, and Mr is the Majorana mass matrix 
for the right handed neutrinos. It is difficult to look for lepton number violation in this case due to 
the singlet nature of the right-handed neutrinos, which does not allow them to be easily produced. 
See Ref. for details. 

• Type II seesaw mechanism j4j: An SU(2) scalar triplet is introduced, A ~ (1, 3, 1), and the neutrino 
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masss matrix is given by M. J = h u v&. Here, h u is the Yukawa coupling between the leptons and 
the triplet, and i?a is the vacuum expectation value of the neutral component of the triplet. In this 
scenario one could find spectacular signals at the LHC from the decays of singly and doubly charged 
components of the Higgs triplet into leptons, i.e. H ± — > e^v and ff ±=t — > efef ID- 

• Type III seesaw mechanism O : It is also possible to generate neutrino masses at tree level by inuo- 
ducing at least two extra fermions in the adjoint representation of SU (2), p ~ (1, 3, 0). The mass 
matrix for neutrinos is similar to the Type I case, where one replaces Mr by M p , the Majorana mass 
matrix for the fermionic triplets. For the testability of the Type III seesaw at the LHC see Q. 

These are the simplest mechanisms for generating neutrino masses at tree level since they involve only the 
addition of one new representation to the minimal Standard Model. It is important to mention that all these 
scenarios can be realized in the context of grand unified theories (GUT's) based on SU(5) and SO(W) 
gauge symmetries. The realization of the Type III seesaw mechanism in the context of GUT's is special 
since it always leads to a hybrid scenario: Type I plus Type III seesaw (SO [3. See Ref. ifTOll for a review 
of the different seesaw mechanisms. 

The simplest mechanisms for the generation of neutrino masses at one-loop level are 

• Zee-Mechanism ifTTTl : Neutrino masses can be generated at one-loop level with the addition of two 
extra Higgs bosons: a Higgs singlet 5 ~ (1, 1, 1) and a Higgs doublet H2 ~ (1, 2, 1/2). In this 
scenario it is difficult to look for lepton number violation because one has only singly charged Higgs 
bosons which do not allow for the like-sign dilepton signals. 

• Colored Seesaw [1]: In this mechanism the "seesaw" fields live in the adjoint representation of 
SU(3). Adding one type of fermionic representation, p ~ (8, 1,0), and one colored scalar, S ~ 
(8, 2, 1/2), it is possible to generate neutrino masses at the one-loop level, as shown in Figure [I] If 
the fermionic fields are light, with mass below or around the TeV scale, one can hope to search for 
lepton number violation through the production of colored states which can be easily produced via 
the strong interactions at the LHC. This is a novel and new possibility which we will explore in this 
paper. 

It is important to mention that in the context of supersymmetric theories with R-parity violation one has 
many interesting predictions for lepton number violating decays which can be found at the LHC. For a 
review see Ref. lfT2l . 
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FIG. 1: Feynman diagram for neutrino mass generation in the colored seesaw mechanism [1 J introducing a scalar 
colored octet S ~ (8, 2, 1/2) and two fermionic octets p a ~ (8, 1, 0). Xhs is the quartic scalar-Higgs coupling and 
Y v the neutrino Yukawa couplings. 

III. COLORED SEESAW MECHANISM 

The Colored Seesaw Mechanism HI introduces two new types of fields: a scalar color octet S ~ 
(8,2,1/2) and a fermionic octet p ~ (8,1,0). Neutrino masses are generated at the one-loop level as 
shown in Fig.[TJ It is important to emphasize that this model is the simplest one of this type at the one-loop 
level which does not require an extra symmetry. The scalar octet has the same quantum numbers, aside 
from color, as the SM Higgs boson therefore leading to a quartic interactions between the two fields and 
couplings to the Standard Model quark fields |[T3l ; the p has the same quantum numbers as the gluino in 
supersymmetric models. 

The relevant interactions (here we write just one possible quartic interaction for simplicity) needed to 
generate neutrino masses are given by, 

-£ u = Y u l T Cia 2 S p + ^ M p Tr (p T C p) + Xhs Tr + h.c. , (1) 

where Y v are the lepton doublet Yukawa couplings to the seesaw fields, Xhs is the quartic scalar-Higgs 
coupling, and 02 is a Pauli matrix and C the charge conjugation operator. Integrating out the seesaw fields 
leads to the usual dimension five operator for neutrino masses. All fermions have left chirality, l T = (u, e), 
S T = (S + , S°), and the trace is over color indices. In order to generate neutrino masses in agreement with 
experiments at least two copies of S or p are needed. Notice that this scenario is consistent cosmologically 
since the scalar octet can decay into SM quarks (T). Finally, we would like to point out that this mechanism 
can be realized in the context of the Adjoint SU (5) |9] grand unified theory since p C 24 and S C 45//. 

In the following we focus on the simplest model with two copies of the new fermions and one scalar 
colored octet. Working in the mass eigenstate basis for the two new fermions p a , a = 1,2, the 3x3 
neutrino mass matrix reads as, 

M ij = yia yja ^HS_ y2 j ^ ^ ^ (2) 
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where the loop integration factor, / (m Pa , ms), is given by 



I{m Pa ,ms) = rn 



(m| - m 2 pa + m 2 p Jn(m 2 pa /m 2 s )) 



Pa 



\m Q — m 



(3) 



l S ""Pa 

With this minimal number of new fields there are two massive and one massless neutrino. If we stick to the 
maximal suppresion coming from the loop factor, the scalar mass should be larger than the fermionic mass. 
This scenario provides the best chance to observe lepton number violation since the Yukawa coupling could 
be large. In the limit ms 3> m p the loop factor simplifies and the neutrino mass matrix becomes, 



KAij = Y ia vi a ^ HS „ 2 



m 



2 • 



(4) 



For illustration purposes let us consider m p = 200 GeV and ms = 2 TeV. With v = 246 GeV, we find that in 
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order to get the neutrino mass "scale", ~ 1 eV, the combination of the couplings has to be Y 2 \hs ~ 10 
Notice that the simultaneous presence of the Yukawa term proportional to Y v and the quartic interaction 
proportional to \hs m Eq. ([T]) violate lepton number. Lepton flavor conservation is also violated, even 
when Xhs is very small, by the Yukawa couplings Y v through processes such as [i — > e + 7 lfi"4l[l~5l . We 
will comment on this in Section ITlIDI 

A. Properties of the Scalar Octet 

The extra Yukawa interactions due to the presence of the scalar octet S are given by 



C Y = d R T D 5 f Q L + u R T v Ql S y e xy + h.c, 



(5) 



where x, y are SU (2) indices and 



S 




ArpA 



ST 



with A = 1, ... ,8 and T A being the SU (3) generators. In the physical basis the Lagrangian reads, 
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Cv =d 



P L [D R T D U L ) - P R (d{t\jUr 



S u + u 



Pr[uIt ] d Dr) -P L (u R TuD L 



S + d 



— i—j=d 
V2 



P L ( D R T D D L ) + P R [D{t ] d Dr 



P L (D R T D D L ) - P R [D\f n D R 



d+ ^ku 

a + i—=u 
V2 



P L [ U R T V U L ) + P R (uItIjUr 



P L [ U R TuU L ) - Pr (ulv\jU R 



u 



(7) 
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Here u and d are respectively the SM up- and down-type quark fields and Ul, U r , Dl and D R are the 
matrices that diagonalize the quark mass matrices. denotes the charged octet scalar and S R 1 are respec- 
tively the CP-even and CP-odd neutral scalars. Assuming minimal flavor violation [13], which we do for 
the rest of this work, means that 

IV = vuYu and T D = VD Y D . (8) 
In this case the physical interactions are 

,MFV t/e. _ _ t/T D_„„T/-t 



C Y =—d\^P L r] D m D V^ KM -P R r] U V ( l KM m u jS u 

H u (Prt] D V C KMin D - P L rju mjjVam) S + d ®) 

v 

+ rjn S R dd + m b B uu + irju Srd^d - irju bjwy 5 u, 

v v v v 

in terms of the quark masses my, mr> and the Cabibbo-Kobayashi-Maskawa matrix Vckm- i]u an d f]D are 
parameters that describe the strength of the scalar couplings to matter. For a list of the Feynman rules see 



Fig. 13 in Appendix [A] 

The masses of the charged state and the neutral members of the scalar octet depend on the details 
of the potential describing their self-interactions. The vacuum expectation value of the SM Higgs boson 
causes a mass splitting between the octet scalars |[T3l . If the splitting is sufficiently large, decays such 
as — > S^W^ are allowed, where S° denotes one of the neutral scalars. However, the mass squared 
splitting is on the order of the electroweak scale squared and will only lead to significant mass splittings 
for light scalars, which are disfavored due to their contribution to b — > sj. Therefore, it is likely that the 
charged scalar will decay predominantly via S + — > tb due to the couplings in Eq. ^ and we adopt this 
assumption for the rest of this study. See Ref. lTT6l for the study of the production mechanisms at the LHC 
and the properties of their decays. 

B. Properties of the Fermionic Octets 

In order to derive the Feynman rules for the fermionic octet we rewrite Eq. ([T]) in a more traditional 
form: 

- C v D YT l u ^ 2 Tr (St p Ra ) + \ Mp Tr (p T Ra C p Ra ) + h.c, (10) 

where the fermionic mass matrix M p can be assumed to be diagonal without loss of generality. Now, 
defining a Majorana fermion as 

f^\ (pr)C \, (11) 

PR 
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the four-component Lagrangian with the mass eigenstates of S has the form: 



- ^= (^MNS^^^P^Tr^ ^) + -J= (rjFp* MN s) a .Tr(^^a) ^ 

- yj Q P fi Tr (5" F Q ) - C* Tr (5+ F a ) P L l~ + Tr F Q F a . (12) 



The corresponding Feynman rules are given in Fig. 13 in Appendix \K 



We note that collider experiments allow to set a conservative lower bound on the mass of the fermionic 
octets of about 200 GeV, as for the gluino in supersymmetric theories. 

C. Constraints from Neutrino Physics 

Since the properties of the fermionic octets depend on the neutrino Yukawa couplings and thus on the 
neutrino mass hierarchy, it is helpful to express the relevant Yukawa couplings as a function of the leptonic 
mixing angles and neutrino masses. Starting from Eq. ([2]), we write the neutrino mass matrix M. v as 

M u = Y u (M ES )- 1 Yj v 2 , (13) 

where M ES is a 2 x 2 diagonal matrix, 

m es = 16*1 diag (/ ( m , ms) , I {m^ms))- 1 ■ (14) 



and I (m Pa , ms) is the loop integration function expressed in Eq. ([3]). The three light neutrino masses can 
be expressed in the following way, 

m u = Vp MNS M u Vp MNS , (15) 

where m v = diag(mi, iri2, m^) and Vpmns can be taken as the leptonic mixing matrix for the three gener- 
ation of light neutrinos without the loss of generality. Working in the basis where the charged lepton mass 
matrix is diagonal, one finds for the neutrino Yukawa couplings Y v , 

Y v = l Vpmns my 2 Q (M Eff )V2. (16) 
v 

Here we are using the Casas-Ibarra parametrization, where Q, is a complex matrix satisfying the orthogo- 
nality condition Q T £l = 1 iflTl . 

In order to understand the constraints coming from neutrino physics let us discuss the relation between 
the neutrino masses and mixing. The leptonic mixing matrix is given by 



^PMNS 



' C12C13 C13S12 e" 1<5 si 3 \ 

-Cl2Sl3S23e l6 - C23S12 C12C23 ~ e 1<5 Sl2Sl3S23 C13S23 
\ S12S23 - e 1<5 Ci 2 C 2 3Sl3 -C23Sl2Sl3e 1<5 - Ci 2 S 2 3 C13C23 / 



xdiag(l,e i4, / 2 ,l) (17) 
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FIG. 2: Feynman diagram for the contribution to \i — > e~f in the colored seesaw model. 



where Sjj = sin dij, Cij = cos 6^ , < Oij < tt/2 and < 5, < 2ir. The phase 5 is the Dirac CP phase, 
and is the Majorana phase. The experimental constraints on the neutrino masses and mixing parameters, 
at 2a level GU, are 

7.25xlO" 5 eV 2 <Am^ < 8.11 x 1(T 5 eV 2 , 
2.18xl(T 3 eV 2 < lAmfj < 2.64 x 1(T 3 eV 2 , 

0.27 < sin 2 #12 < 0.35, (18) 

0.39 < sin 2 6*23 < 0.63, 
sin 2 6>i3 < 0.040, 

and J2i m i < 1-2 eV IT91 . Following the conventions, we denote the case Am^ > as the normal 
neutrino mass hierarchy and Am^ < the inverted hierarchy, i. e.we have 

m-i = 0, 7772 = \j Am^, 7773 = yj AWjj | (NH), 

7772 - \/ .Allrr, | -|- ^777^ , 7773 = (IH). 



(19) 



III ] — \/ \Am^, 7772 = \J I Amj] I + A?772i, 777-3 = 

The 17 matrix takes the well-known form corresponding to the Type I seesaw case with two right-handed 
neutrinos [20]. Expressed in terms of an angle oj, it reads 

/ n n \ / y/\ ->.V \ 



n 



NH 



-OJ 



V 



OJ 



OJ' 



J 



-OJ 



\ 



OJ 









(20) 



/ 



in the normal (NH) and inverted (IH) hierarchy, respectively. Here we focus on the range, — 1 < oj < 1, for 
simplicity. 



D. Constraints from Lepton Flavour Violation: /j -> £7 

Neutrino oscillations indicate the presence of lepton flavor violating operators which lead to various rare 
decays, the most stringently constrained being fx — > ej. The colored seesaw has a contribution which can 
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potentially disagree with experiments (see Fig. [2]). The branching ratio for this process is given by 



BR (p ej) 



r(/x -> ej) 
T(fi — > evv) 



3(Xem 



(21) 



where x a = m 2 Pa / m 2 s and 



l-6x + 3x 2 + 2x 3 -6x 2 ln2; 



(22) 



12 (a; -If 

Here we have used T(fi — > evv) = G 2 7 m 2 l /(1927r 3 ). Notice that our result agrees with [15] up to the 
normalization for the Yukawa couplings, but disagrees with |[T4l . It is clear from the above that this branch- 
ing ratio is very sensitive to the Yukawa couplings: BR (/j, — > ej) ~ (^) 4 - Only for order one Yukawa 
couplings this process does conflict with the experimental bound of BR(/i — > ej) < 1.2 x 10 -11 ED , 
which will be improved by two or three orders of magnitude in future experiments. This corresponds to 
^HS ~ 10~ 10 for the quartic coupling between the octet scalar and the SM Higgs. While one can claim that 
this is inline with the spirit of the seesaw mechanism (order one Yukawa couplings) there is no argument 
for why the quartic coupling should be so small, although such a small value is protected by global U{1)l 
lepton number. Nevertheless, the interplay of the low energy constraints with the collider signatures can be 
exciting in this part of the parameter space. In Fig. [5] we show the predictions for BR(^ — > ej) versus the 
values of sin #13, for 1 TeV < = mF 2 , fns < 5 TeV in the case of NH and IH, respectively. One finds 
that the quartic coupling is restricted to 



\hs > 10" 



(23) 



in order to comply with the experimental bounds on the rare decay. In the case of degenerate fermions, 



Eq. (21 1 simplifies (the u dependence drops out). Using Eq. ( 16 1 to express the Yukawa couplings in terms 
of the neutrino parameters and the quartic coupling Xhs ( see Appendix [B] for Yukawas explicitly given in 
terms of neutrino parameter) and expanding in S13 shows that the rate in the NH goes as 



BR (fx — > ej) oc 



C12 C23 S12 



\m\ x + si 3 s 2 3e 



\ 2 ' 



(24) 



This a sum of two terms which are both suppressed: the first by the solar mass parameter and the second 



by S13. This leads to a cancellation at around S13 ~ 0.09 and 5 — * = 
left-hand panel of Fig. [3] In the case of IH the branching ratio goes as 



7T, corresponding to the dip in the 



BR (/x — > ej) oc 



-S13 S23 e 



~V 1 

m 31 + o C 12 s 12 C23 



Am 2 21 



A HS 



(25) 



V Am ii 

where a further expansion in the solar mass scale has been conducted. Once more we have the sum of two 
suppressed terms with a zero at around S13 ~ 0.007 and 5 = 0, again reflected in the right panel of Fig. [3] 
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FIG. 3: BR(/i — > ej) as a function of sin #13 for 1 TeV< ms,mp <5 TeV with mp 1 = mp 2 . In the left panel 
normal hierarchy (NH) and in the right panel inverted hierarchy (IH) is considered. 

Since the neutrino Yukawa couplings define the decay length of the fermionic octets, as we will show in the 
next section, the above results thus define a lower bound on the decay length. 



IV. COLORED SEESAW AT THE LHC 



In this article we focus on the pair production of the fermionic octets to understand the predictions for 
the lepton number violating channels with multileptons. The production mechanisms for scalar octets do 
not lead to lepton number violating signals and have been investigated in previous studies lfT6l . 



A. Decay Modes of the Fermionic Octets 

Lepton number violation in the context of the colored seesaw model originates from the decays of the 
fermionic octets, F a . Depending on the mass configuration of the octets, different decay channels dominate. 
\fmp a > m s ± (and mp a > m s o ), the dominant decay modes are two-body decays into a scalar octet 

R/I 

state and a lepton, 

F a — > < l \ S and F a -»■ J Sr/i . (26) 
S + y Ui Sx/j 

Here denotes a charged lepton and i = 1, 2, 3 is the generation index. However, if m s ± > mp a > TUt, 
the fermionic octets predominantly decay via an off-shell (charged) scalar octet into a lepton and a heavy- 
flavor quark pair, 

(27) 

ertb 
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The three-body decay can also be mediated by a neutral scalar octet state, if m q o > nip > 2m t , 

J R/I 



Vitt 

v. ti 



(28) 



Similar decays into a bottom- or light-flavor quark pair are suppressed by the small quark Yukawa couplings 
and contribute negligibly (given the assumption we have adopted of MFV). 



1. Two-body decays 

The decay width of a general two-body decay — > 1 2 follows directly from the spin and color averaged 
squared matrix element \Ai\ 2 [22], 

T(0^ 12) = \M(0-> 1 2) I 2 , (29) 

87r7TT,Q 1 1 

where p cm is the momentum of the final-state particles in the rest frame of the decaying particle, p t 



cm 



l^r K - {mi + m 2 ) 2 ] [m 2 , - (mi - m 2 ) 2 ] . 

Here we consider the decay of a fermionic octet F a of mass tof q into a charged lepton (neutrino) and a 



charged (neutral) scalar octet state of mass m s ± (m s o ), cf. Eq. (26). We neglect the lepton mass in the 

R/ 1 



decay rate calculations. Using the Feynman rules given in Appendix [A] one finds 

\M(F a -+ lfS±)\ 2 = 1 |yj Q | 2 (ml a -m| ± ), 
M{F a -+ ^5° /7 )| 2 = 1 \{V? Mm Y u ). a \ 2 (m 2 Fa - m^] 
and thus the partial decay rates are 



(30) 



128vr \ p a 

let 

F(F a -> ui S° R/I ) = -L | (V? MNS Y U )J (m\ a - ml 0R ) * 



(3D 



6 ^ 

-r a 

In the following we neglect the mass splitting between the scalar octet states, as favored by b — > sj, and 
assume mj = m s ± m s o m m s o. The total two-body decay rate T2 tot of fermion F a is obtained by 

It I 

summing over all three lepton generations and by including the charge conjugated decays, 

r 2 ,tot(^) = 3^ £ |>T | 2 (m% a - m|) 2 -L-. (32) 

i Fa 

The neutrino mixing matrix Vpmns cancels out in the sum over all three neutrino generations. The resulting 
two-body branching ratios are remarkingly simple and do not depend on the precise mass configuration of 
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the fermion and scalar states, 

BR(F a -> £~s + ) + BR(F a -> #s~) = ^ |y7| 2 / ]T |y; Q | 2 , 

* (33) 

BR(F Q -> = i 
where the sum of the branching ratios to charged leptons is then also one half. 

2. Three-body decays 

The decay width of a general three-body decay — > 1 2 3 can be parameterized as follows l22l . 

r(0 1 2 3) = TT^^r^ / W ^ / 23 d P23 |^(0 ->• 1 2 3)| 2 , (34) 
(27r) d 32m^ J v mm J p gn 



with pjj = (pi + pj) 2 and pi and m, denoting the momentum and the mass of particle i, i = 0, 1, 2, 3. The 
integration limits are given by, 



# = (mi + m 2 



,2 



P?£ x = (m -m 3 ) 2 , (35) 
^n/ max = (£3 + El? _ ^ Ef _ m , ± _ ^ 

where E\ = (pi2 — m\ + m 2 .) / (2 v /pi2) and = (m 2 , — pi2 — m 2 )/ (2^12) are the energies of particle 
2 and 3 in the p\2 rest frame, respectively. 

The squared matrix element for the three-body decay Eq. ( |27| ) reads, in terms of the invariant p23> 

\M(F a -» C (5)| 2 = j |Y»f ,36) 



The three-body decay Eq. (28 ) can be mediated by an off-shell or Sj scalar. The summed matrix element 
is obtained as follows, 

\T7fjp , ++ "M 2 1 \(vT v \ I 2 ^U m t , 2 ^/(P23-4m 2 ) ^23 1 

LM(-F a -)> I/i*t) = — (VpMJVS^) ,•„ "^V" ( m F a ~ P23) < 7 0-^2 + 



16 ir^^^ fa | ^2 V* a W 1 (;p2 3_ m 2 o) 2 (p 23 - m 2 o) 2 



(37) 



For equal masses, 7715 = m s o^^, the matrix element becomes, 



\M(F a -> ^ = I | ( ^ ws n) .J 2 ^ { ( '" L ^ ( :i ) ; 2 '"' ) } • 08, 

Since the top-quark mass cannot be neglected, the matrix elements and the partial widths for the charged 
lepton channel and the neutrino channel have sligthly different kinematics. In contrast to the case of two- 
body decays, the resulting branching ratios for the three-body decays thus do have a residual dependence 
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FIG. 4: Two-body branching ratios in the case of degenerate fermion masses mj?i = m^ 2 as a function of sin# 13 
both for normal hierarchy (NH, left) and inverted hierarchy (IH, right). Full scan over the leptonic mixing parameters 
according to Eq. ( |18) , The bands for the \i- and r-rates overlap, where the band width is mainly given by variations 
in sin #23. 

on the masses of the new color octet states. This mass dependence is small however and for the study of the 
branching ratios we will focus on the flavor patterns as described by the two-body decays. It is important to 
mention that when the mass of the fermionic octets are in the range, m t < mp a < 2m t , the F a can decay 
only into charged leptons and the three-body decays into neutrinos is kinematically forbidden. Therefore, 
in general, the above given results define lower bounds on the rates for decays into charged leptons. 

3. Numerical Results 



We begin by focusing on the simplest case where the fermionic octets are degenerate in mass. In this 
scenario the branching ratios are independent of the unknown matrix. Therefore, the results will only 
depend on the neutrino mixings and phases. In Fig. [4] we show the total branching ratios versus sin #13 
with a scan over all other neutrino parameters in the NH (left panel) and IH (right panel), respectively. 
Remarkably, the gross features of the branching ratios are independent of the neutrino parameters as is 
the value for the branching ratio to electrons. This property can be understood by expressing the Yukawa 
couplings in terms of neutrino parameters (see AppendixlB]) and expanding in S13. For the NH: 



Am 2 31 , , 



T(F^eX) oc \ si 2 ^Am 2 2l + si 3 , 
T (F —?■ fiX) oc (c 2 12 4^ + s 2 2 ^AmlA , 



r (F -> tX) oc c{ 2 s 2 23 JAm 2 21 + c 2 23 J Am 2 31 



(39) 
(40) 
(41) 



14 



From these expressions we see that in the electron branching ratios, the larger atmospheric mass is sup- 
pressed by s\ 3 whereas in the muon and tau channels this term has a factor of about a half (for tri-bimaximal 
mixing) leading to the order of magnitude hierarchy between these branching ratios. This is a clear pre- 
diction of this degenerate fermionic mass scenario. Furthermore, since the leading term in the electron 
branching ratio is the s^Am^, which depends on more accurately measured parameters, the branching 
ratio itself can be predicted more accurately. 
In the IH case one has: 

r (F -> eX) oc y/Aml v (42) 
T (F ^ »X) <x 4 3 ^J Am 2 31 , (43) 
r (F -> tX) oc 4j^/a™L ( 44 ) 

where the solar mass scale has been neglected. Here we see then that all branching ratios are proportional 
to the atmospheric mass scale with a slight suppression for the muon and tau channels, about a half for 
tri-bimaximal mixing. Again the parameters in the electron branching ratio have small uncertainties leading 
to a more certain prediction for its value. The upshot of this is that for the case of degenerate fermions, 
there is a definite prediction for the branching ratios, and measuring these would distinguish between the 
two possible neutrino spectra: 

BR (F -> fiX or tX) » BR (F -> eX) => NH, (45) 
BR (F -> eX) > BR (F ->■ yX or tX) => IH. (46) 

For the more general scenario (non-degenerate Fs) we start by assuming tri-bimaximal mixing in the 
neutrino sector for illustrative purposes, (9^ = 0, sindJ 2 B = l/y/3, sin 9^ = l/y/2). The results are 
shown in Fig. [5] referring to F\ (F2) in the upper (lower) panels. Already we see that the dependence on oj 
is quite important, even in this simplified case for neutrino mixing, since the branching ratios vary widely 
with uj and all of them go to zero at some specific oj value. However, general statements can still be made: 

• In the NH case either the /^-channel or the T-channel always dominates over the e-channel. 

• In the IH the BR(F a — > fiS) and BR(F a — > tS) coincide. The e— channel can be dominant for 
some values of uj. 

Next we scan over all the neutrino parameters in order to understand their effect on the fermionic decays. 
We show the results in Fig. [6] for the branching ratios of F\ and F2 versus oj. This is of course even more 
complicated than Fig. [5] and while we can still state that in the NH either the /^-channel or the r-channel 
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FIG. 5: Two-body branching ratios for mp > ms as function of to both for normal hierarchy (NH, left) and inverted 
hierarchy (IH, right) in the case of non-degenerate fermion masses. Upper (lower) plots refer to Fi (F 2 ) decays. 
Results do not depend on the precise mass configuration. The parameter uj determines the neutrino Yukawa coupling, 



see Eq. ( 20 » and Eq. ( 16 1. Central values for the neutrino mass differences and tribimaximal mixing in the leptonic 
mixing matrix are chosen without CP phases. 



always dominates the e-channel, we can no longer make the statement that in the IH the BK(F a — >• /xS") 
and BR(Fq, — > tS) coincide, due to varying S23 and S13 away from the tri-bimaximal values, which treat 
the muon and tau neutrinos equally. It is important to mention that if this model is realized in nature 
and one finds the color octet fields at the LHC, and the neutrino spectrum is discovered in future neutrino 
experiments, then in principle one could determine the corresponding value of u from the comparison with 
Fig-© 



B. Production Mechanisms and Lepton Number Violation 

In this section, we first present the cross section for the pair production of fermionic octets at the LHC. 
Second we will discuss the decay length. Finally, we combine our findings with the results of the decay 
properties. We focus on the signals for lepton number violation due to two final-state leptons of the same 
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FIG. 6: Same plots as Fig. [5] but now the leptonic mixing parameters are varied between the experimental bounds 



given in Eq. ( 18 1. Nonzero CP phases are included, < <5 < 2ir. For the normal hierarchy, the majorana phase is 



varied between < $ < it only for illustration purposes 
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FIG. 7: Parton-level Feynman diagrams for the pair production of color octet fermions at lowest order, via gluon- 
gluon fusion (left) and quark-anti-quark annihilation (right). 



electric charge. 



1. Production offermionic octets at the LHC 



The fermionic octets in the colored see-saw model can be produced at hadron colliders via the strong 
interaction. The dominant production mechanism proceeds at 0{a 2 s ), via the two partonic channels gluon- 
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gluon fusion and quark-anti-quark annihilation, 

g(pi) g(P2) — > F a (p 3 ) F a (p A ) and q(pi) q(p 2 ) — > F a (p 3 ) F a (p 4 ), (47) 

where q = u,d,c,s,b can be any of the light-flavor quarks. The corresponding Feynman diagrams are 
shown in Fig. [7] Mixed pairs of fermions, F\F%, cannot be produced since the gFF and ggFF vertices are 
diagonal with respect to the fermionic eigenstates. We parameterize the cross sections in terms of the usual 
Mandelstam variables, 



s = (pi+P2) 2 , i = {Pi ~ P?,) 2 , u = (p 1 -p i f 



(48) 



The differential partonic cross sections for the subprocesses are obtained from the spin- and color-averaged 
squared matrix elements, 



I 2 

dOgg,qq{s) = ^~^2 \Mgg,qq(s, t, u) \ dt, 



(49) 



which can be written as follows, 



\M gg (s,i,u)\ 2 = 187r 2 a^ 



1 - p(t - m F J(u - m Fa ) 



2 + 



4ml s 

fa 



i-r-7 ,„ 2 ,s,2 6An 2 a 2 s 1 

\Mqq{s,t,U)\ = ^ 

o S 



2m 2 F J + (t - m 2 Fa ) 2 +(u- m 2 Fa f 



s 



(t - m 2 F J{u - m 2 F J 



(50) 



Here, a factor 1/2 has been taken into account because of the identical particles in the final states. Note 
that the results for pair production of the fermionic states are the same as those for gluino pair production 
in supersymmetry with decoupled squarks 11231 . 

At the hadronic level, the cross sections are obtained from the partonic ones by the convolution, 



dapp^ FaF (s) 



1 f dC pp dC pp 



TO 



dr * dr f' 



da gg (s) + 



(51) 



with r = s/s, s being the hadronic center-of-mass energy squared, and tq = Am 2 Fa /s is the production 
threshold. The parton luminosities are given by, 



aL -ab 

dr 



i r L dx r 

1 + $ab Jr x 



(52) 



where the parton distribution functions (PDFs) f a /A( x , I 1 ) parameterize the probability of finding a parton 
a inside a hadron A with faction x of the hadron momentum at a factorization scale p. 
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FIG. 8: Hadronic cross section for the pair production of color octet fermions within the framework of the colored 
seesaw model at the LHC with ^/s = 7TeV (solid) or s/s = 14TeV (dashed), as a function of the mass of the 
produced fermion, mp a . The cross sections for the dominant background, tttt, are also shown, taken from |26|. 

The numerical results are shown in Fig. [8] both for ^/s = 7TeV and ^/s = 14TeV, as a function of the 
mass of the produced fermionic states. The PDF set CTEQ 6L G5l is chosen and the factorization scale is 
set to \i = niF a - Only for small momentum fractions x and thus for light fermion masses, the qq channel 
is suppressed from the PDFs and less important than the gg channel. At a 14 TeV collider, the fraction 
of qq and gg initial states is about 25% and 75%, respectively, for a fermion of mass m,F a = 1 TeV. For 
yfs = 7 TeV the qq channel dominates over the gg channel for mp a > 1.2 TeV. The full hadronic cross 
section falls steeply with increasing fermion masses, from a « 300 pb for m_p a = 200 GeV to a « 1 fb for 
mp a = 1 TeV in the early stage of the LHC at y/s = 7 TeV. For illustration we also show in Fig.[8]the cross 
section for the leading background, tttt, as we will discuss below. 

2. Decay length and stable color octet fermions 

The possible decays of the fermions depends on the ratio of the color octet masses, the neutrino Yukawa 
couplings, and the parameter Xhs- In Fig. |9]we show their decay lengths, scanning over — 1 < oj < 1, 
for ms = 1 TeV in case of normal hierarchy (NH, left in Figj9]) and inverted hierarchy (IH, right in Fig. 
[9]) versus the fermion mass. Since we sum over all lepton generations, the dependence on the neutrino 
parameters disappears due to the unitarity of Vpmns> hence there is no need to scan over these parameters. 
This can be seen explicitly by expressing the Yukawas in terms of the neutrino parameters and u and 
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FIG. 9: Decay length of the fermion for mg — 1 TeV in case of normal hierarchy (NH, left) and inverted hierarchy 
(IH, right). For mp < ms the fermion decays via three-body decays, for mp > ms the fermion undergoes two-body 
decays. We scan over the parameter — 1 < to < 1. The decays for F\ and F^ only differ for specific values of u and 
agree when scanned over the full range of ui. 



expanding in S13 (see Appendix [B]). The decay length in the NH and IH is then proportional to: 

-1 



d NH oc 
dm oc 



Am 31 



-1 



(53) 



respectively. In the inverted hierarchy case it was possible to make a further expansion in the solar mass 
scale which led to the cancellation of the u dependence at this order and is clearly reflected in Fig. [9] 
Although the u dependence is not strong in either spectrum. For mp < ms the fermion decays via three- 
body decays, for mp > ms the fermion undergoes two-body decays. It is clear from the plot that the 
parameter space allows for everything from prompt decays to displaced vertices and even decays outside 
the detector (large Xjjg small fermion mass). 

In the latter case where the fermions decay outside the detector, they exit the detector as so called R- 
hadrons, i.e. they hadronize. Such scenarios have been studied before, especially in supersymmetry where 
the gluino can be (meta-)stable [24]. This scenario is analogous to ours since the gluino's quantum number 
are the same as our F's. While such long lifetimes are discouraging for the observation of lepton number 
violation at colliders, the associated signals are still a spectacular indicator of new physics. 



3. Signals for lepton number violation 

The most important production mechanism for lepton number violation is the QCD pair production of 
fermionic octets. Each of these can subsequently decay into a charged lepton and either an on-shell or off- 
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shell charged scalar octet which further decays into a top-bottom pair. In order to preserve the signature of 
same-sign dileptons, we must further require the top to decay purely hadronically into three jets, t — > 3 j. 
Then the signals of interest are CQ 

PP -> F a F a — > e ± e ± 8j; 

— ► f f 8 j; (54) 
— ► e ± 8 j. 

The tau channel is also very important to be help distinguish between the different neutrino spectra. The 
event reconstruction will be much harder to carry out at the LHC. We continue by focusing on the simple 
case where the fermions are nearly degenerate therefore allowing more predictivity. The number of events 
N per 10 fb -1 for each of the signals is then given by 

N(£f£f) = Q{m Fa - m s ) N 2hoAy (£±£f) + Q(m s - m Fa ) N 3hody (£f£f), (55) 

with 

N 2h ody(tftf) = 2 °PP->F a F a x BR(F a -> £fS^) x BR(F a -> £fS^) 
a=l,2 

x [BR(5 T -»■ 46)] 2 x [BR(t -4 3 j)] 2 x 10 (56) 
« 2 ^ (Tpp—^ p a p a x BR(F Q -> £f 5 T ) x BR(F a -> £±S^) x 4.6, 

a=l,2 

and 

^3body(^^) = 2 crpp^ Fct x BR(F a £ftb) x BR(F Q -> *±tb) x [BR(t -4 3j)] 2 x 10 

a=l,2 

« 2 ^ a PP ^ FaFa x BR(i ? a -> *fib) x BR(F a -> *±tb) x 4.6, 

a=l,2 

(57) 

where the first ©-function indicates two-body decays and the second three -body decays, the factor of 2 
comes from the sum over the charge conjugated decays, we have used BR(i — y 3 j) ~ 0.68 and the cross 
section is in fb. As an order of magnitude estimate, for 500 GeV fermions and a branching ratio to muons 
at about 20% at the LHC for 7 TeV center of mass energy one would expect about 180 like-sign dimuon 
events, before accounting for signal efficiency, which could play an important role. Requiring two iso- 
lated like-sign leptons and eight jets, the leading background comes from tttt production where one pair of 
like-sign tops decay semi-leptonically and the rest decay hadronically yielding the signal 1^8 j +missing 
energy. The cross section for tttt is quite small, about 2 (11) fb at y/s = 7 (14) TeV [26]. Addtitional back- 
ground contributions arise from three tops plus jets, which is subleading due to the electroweak production 
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FIG. 10: Rates for like-sign dilepton events per 10 fb~ from the pair production of color octet fermions at the LHC 
with ,/s = 7 TeV in the case of normal neutrino hierarchy (NH, left) and inverted hierachy (IH, right). We scan over 
the leptonic mixing and lu = (upper plots) or — 1 < lu < 1 (lower plots). The fermions are assumed to be nearly 
degenerate in mass, mp 1 sa rap 2 . The common scalar mass is set to ms — 2000 GeV for illustration purposes, but 
the results are not very sensitive on the scalar mass. 



mechanism (0.5 fb at 14 TeV II26II ). ttbb in comparison has a larger inclusive cross section (about 2600 fb at 
14 TeV E71 ). but the signal rate depends strongly on the fake rate for an isolated lepton resulting from a b 
decay. Also, if we stick to the possibility that the top quarks can be reconstructed, see for example Ref. [28 ], 
and the simultaneous presence of two same-sign leptons, ttbb channel is not a relevant background in our 
study. Of course a more detailed study must be conducted to understand the observability. 

and 



We present the signal and background in Figs. 10 



11 



at the LHC with 



7 and 14 TeV, respec- 



tively. In both panels, the bands are due to the scan over the neutrino parameters. The upper panels are with 
the choice of u = 0, and in the lower panels we scan over — 1 < uj < 1 to indicate the possible range of 
the signal rate. The common scalar mass is set to ms = 2000 GeV, i. e. the fermions decay via three-body 
decays. As explained above, the results for parameter regions with two-body decays look very similar. We 
see that the number of signal events is larger than the SM background up to fermion masses of about 1 TeV 
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FIG. 1 1 : Same as Figure 10 but for the LHC with Vs = 14 TeV. 



for the dimuon channel at this level of analysis. 

The rates can be understood by studying the results of Fig. [5] (and [6]>. For the specific case of u = 0, 
for the normal neutrino hierarchy (left), the electron and the muon rate of F\ decays are identical while 
the electron channels for F2 decays are zero. As a consequence, only the pair production of F\F\ can lead 
to e ± e ± and e^p^ signals, and they have the same rate, while both F\F\ and F2F2 contribute to \x 
signals. For the case of inverted hierarchy (right) and uj = 0, F2 production contributes equally to e and /i 
rates and F\ predominantly to e signals. 

In both the top panels and the bottom panels, the neutrino spectrum has left its mark on the ratios 
of the e e and n \i number of events. Specifically in the normal hierarchy the number of muon events is 
larger than the number electrons events and vice versa in the inverted hierarchy. Then we have verified the 
earlier statement that in addition to seeing the like-sign dileptons, in the special case where the fermions 
are degenerate comparing the number of muon events to the number of electron events can also reveal the 
neutrino spectrum. 

So far, we have not included any experimental acceptance. The additional cuts will undoubtedly reduce 
the signal rate by a significant factor. On the other hand, the kinematical requirements in the events, such 
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as no large missing energies, mass peaks at mp a « m(ltb) and possibly at ms ± ~ m(tb), would help the 
signal identification and further background suppression as well. We leave the detailed simulations to future 
experimental studies. For the study of the trilepton channels in see-saw models see Ref. [29]. 



V. SUMMARY 



There is an interesting possibility that the fields generating neutrino masses at one-loop level are scalar 
and Majorana fermionic color-octets of SU(3)c, the colored seesaw mechanism. We reiterated the main 
features of this model, in comparison with the other common seesaw theories. We presented the bounds on 
certain model-parameters from the low energy decay \i — > e^y and from the observed neutrino mass spectrum 
and oscillation mixing parameters. We found that these parameters may have interesting correlations to the 
decay branching fractions of the colored states. 

In this context, we investigated the lepton number violating signals at the Large Hadron Collider from 
the production and decay of the fermionic color-octet states. Due to the QCD strong interaction, these 
states may be produced at the LHC with a favorable rate. We focused on the conclusive lepton number 
violating channel: same-sign dileptons (muons and electrons) channels. We found that for fermionic octets 
with mass up to about 1 TeV, the number of same-sign dilepton plus multiple jets events is larger than the 
SM background, indicating that this might be a very promising signal at the LHC which deserves a more 
detailed study. One may be able to use the fermionic octet decays of the lepton flavor combinations to 
distinguish between the neutrino spectra. This becomes especially predictive in the simplest case when the 
two fermionic color-octets are degenerate in mass. Our results demonstrate the potential of the LHC in 
searching for new physics beyond the SM in connection with neutrino mass generation. 



Appendix A: Feynman Rules 



The Feynman rules for the Yukawa interactions of the fermionic octet are collected in Figure 12 and 



those for the scalar octet in Figure 13 Here a = 1..2 indicates the different generations of F, i = 1..3 the 
generations of the leptons, and upper (lower) case roman letters denotes indices of the adjoint (fundamental) 
representation of SU(3). Vtb is the corresponding entry of the quark mixing matrix Vqkm- 



Appendix B: Yukawa Couplings 



As indicated in Eq. ( [To} the Yukawa couplings can be expressed in terms of the neutrino parameters and 
the matrix Q. This can be illuminating for calculating various branching ratios and we do so explicitly here 
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FIG. 12: Feynman rules for the Yukawa interactions of the fermionic octet. The factor of one half results from the 
trace of the SU(3) generators. 




FIG. 13: Feynman rules for the scalar octet. 

where we have expanded in S13. We furthermore use Eq. [l4]to express the results in terms of Xhs- In the 
following, I a = I(m Pa ,ms) where a = 1..2 (the dimensionality of I a is inverse mass). 
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The partial widths (and branching ratios) for the decays of the octet fermions are proportional to Yukawa 
couplings squared. Since the Yukawa couplings are the sums of two terms with different dependence on oj, 
they always cancel at some value of oj for a zero phase. The total width involves a sum over all three lepton 
generations and as such is independent of the Vpmns parameters: 
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Again, this is true for two or three body decays. The two terms in the electron branching ratio are both 
suppressed, one by the solar mass scale and the other by the sin of #13 so that the other two branching ratios 
dominate due to the atmospheric mass term. These two branching ratios are equal in the tri-bimaximal case 
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which goes to zero at s\s ~ 0.09 and S — * = 7r. 
Inverted Hierarchy 
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The total width in the inverted hierarchy is 
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so that varying oj does not much influence on the width. 
For degenerate octet fermions masses and $7 6 1Z: 
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Here there is a very mild hierarchy between the electron branching ratio and the muon and tau do to S23 ~ 
C23 ~ 0.5 and this is the only parameter that plays an important role. Again the muon and tau branching 
ratios are equal in the tri-bimaximal case. 

The branching ratio for \i — > ej is a bit tricky in the inverted hierarchy since it vanishes in the zeroth 
order expansions used above. Therefore, one must use the full PMNS matrix: 

2 

-■si3S 2 3e -\ii\mti + -C12S12C23 , " • (B8) 
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This goes to zero at around sis ~ 0.007 and 5 = 0. 
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